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Abstract. I review our recent results concerning the metal abundances 
in quasar absorption-line spectra obtained with the UVES/VLT, HIRES/ 
KECK, and STIS/HST spectrographs. The analysis is based on the 
Monte Carlo inversion procedure aimed at the restoring the physical pa¬ 
rameters and the ID distributions of the velocity and density fields in 
quasar absorbers. A functional dependence of the line-of-sight velocity 
dispersion, Uv, on the absorber linear size, L, was found: the majority of 
the analyzed systems follow the scaling relation ~ (Ah L)^'^ (with Ah 
being the total gas column density). This may indicate that the metal 
absorbers are associated with virialized systems like galaxies or their ha¬ 
los. A metal content enhancement (up to Z ~ was observed in 

small-size absorbers (L ^ 0.4 kpc), whereas decreasing metal abundances 
were found in systems with increasing L. The lowest metallicity ([C/H] 
~ —3.0) was detected in the Zabs = 2.917 Lymna-limit system (LLS) with 
A(H i) = 3.2 X 10^^ cm“^ and L ~ 140 kpc. The relative abundance ratio 
of [Si/C] ~ 0.35, measured in this LLS, seems to indicate that the initial 
mass function for the stellar population, produced the observed metallic¬ 
ity pattern, was constrained to intermediate masses, Mup ^ 25Mq. 


1. Introduction 

Metal absorption-line systems in quasar spectra are traditionally used to study 
different physical processes in the intergalactic matter (IGM) over a wide range 
of redshifts, back to the time when the Universe was less than 7% of its present 
age. It is supposed that at z 10 the IGM was pre-enriched by the first 
Population III (Pop III) stars formed from gas with zero metallicity (e.g., Bromm 
et al. 2001; Nakamura & Umemura 2001) and enriched later due to disruption 
of low-mass protogalaxies (e.g., Madau et al. 2001) at 6 ^ z ^ 10, and due 
to ejection of metals from massive galaxies at z 6 (e.g., Aguirre et al. 2001; 
Scannapieco et al. 2002). 

The simulations of collapse and fragmentation of primordial gas clouds sug¬ 
gest that the first generation stars may have been very massive, M* IOOMq 
(Abel et al. 2000; Bromm et al. 2001). However, the calculations by Naka¬ 
mura &: Umemura (2001) produced a bimodal initial mass function (IMF) for 
Pop HI stars with a second peak at 1 — 2Mq. According to these calculations 
the first generation supernova (IOMq < M* < 35Mq) and pair-instability su¬ 
pernova (140Mq < M* < 260Mq) can produce metallicity in the range 10“^ 
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to 10~^Zq. Since the element yields and production rate differ significantly for 
the massive (10 — 35Mq) and very massive (140 — 260Mq) stars (Woosley &: 
Weaver 1995; Umeda & Nomoto 2001; Heger &: Woosley 2002), the Pop III IMF 
can be constrained by measuring the abundance ratios in low metallicity cosmic 
objects. 

Unfortunately, the computational methods usually applied to high resolu¬ 
tion spectra of QSOs lie quite often behind the quality of observational data 
and fail to extract from them all encoded information. The common process¬ 
ing method consists of the deconvolution of complex absorption profiles into 
an arbitrary number of separate components (assuming a constant gas density 
within each of them), which are then fitted to Voigt profiles. However, in many 
cases this procedure may not correspond to real physical conditions; observed 
complexity and non-similarity of the profile shapes of different ions indicate that 
these systems are in general absorbers with highly fluctuating density and ve¬ 
locity fields tightly correlated with each other. The need for more sophisticated 
procedures of data analysis is therefore obvious. 

In recent years, it has been shown that accounting for the correlations in 
the velocity field {mesoturbulence) may change the interpretation of the absorp¬ 
tion spectra substantially (e.g., Levshakov et al. 2000a; Levshakov et al. 2002, 
hereafter LACM). Our first inversion codes, — the Reverse Monte Carlo (Lev¬ 
shakov et al. 1999a), and the Entropy-Regularized Minimization (Levshakov et 
al. 1999b), — have been developed for a model of a stochastic velocity field 
neglecting any density fluctuations. They have been applied to the analysis of 
the H I and D i lines and/or to the metal absorption lines with similar profiles 
when different species trace the same volume elements independently of the den¬ 
sity fluctuations. Later on, we extended this study to the inverse problem for a 
model of compressible turbulence when one observes non-similar profiles of dif¬ 
ferent atoms and/or ions from the same absorption-line system. The developed 
procedure, called the Monte Carlo inversion (MCI), is described in Levshakov 
et al. (2000b, hereafter LAK), and its advanced version in LACM. 

The MCI code has been applied to the analysis of metal absorbers with 
4 X 10^^ cm“^ < A^(H i) < 5 x 10^^ cm“^ from the following QSO spectra; 
J2233-606 (LACM), Q0347--3819 and APM BR J0307-4945 (Levshakov et al. 
2003a), HE 0940-1050 (Levshakov et al. 2003b), HE 0515-4414 (Agafonova et 
al. 2003a), and PKS 0528-250 (Agafonova et al. 2003b). A brief description of 
the obtained results is given below. 


2. Model assumptions and the MCI procedure 

Numerous hydrodynamical calculations performed in the previous decade have 
shown that the QSO absorption lines arise more likely in the smoothly fluctuat¬ 
ing intergalactic medium in a network of sheets, filaments, and halos (e.g., Cen 
et al. 1994; Miralda-Escude et al. 1996; Theuns et al. 1998). A very important 
characteristic of the continuous absorbing medium is that the contribution to 
any point within the line profile comes not only from a single separate area (a 
‘cloud’) but from all volume elements (‘clouds’) distributed along the sightline 
within the absorbing region and having the same radial velocity (see, for details, 
§ 2.2 in LAK). 
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The MCI procedure is based on the assumption that all lines observed in 
a metal system arise in a continuous absorbing gas slab of a thickness L with a 
fluctuating gas density and a random velocity field. We also assume that within 
the absorber the metal abundances are constant, the gas is optically thin for the 
ionizing UV radiation, and the gas is in the thermal and ionization equilibrium. 
The last assumption means that the fractional ionizations of different ions are 
determined exclusively by the gas density and vary from point to point along the 
sightline. These fractional ionization variations are just the cause of the observed 
diversity of profile shapes between ions of low- and high-ionization stages. 

It is well known that the measured metallicities depend in a crucial way 
on the adopted background ionizing spectrum. We started in all cases with 
the Haardt-Madau (HM) background ionizing spectra (Haardt & Madau 1996) 
computing the fractional ionizations and the kinetic temperatures with the pho¬ 
toionization code CLOUDY (Ferland 1997). If the fitting with the HM spectrum 
was impossible, we used another spectra. 

The MCI procedure itself is implemented in the following way. Within the 
absorbing region the radial velocity v(x) and the total hydrogen density nH(a;) 
along the line of sight are considered as two random fields which are represented 
by their sampled values at equally spaced intervals Ax, i.e. by the vectors 
{ui,...,Ufc} and {ni,...,nfc} with k large enough (~ 150 — 200) to describe 
the narrowest components of the complex spectral lines. The radial velocity 
is assumed to be normally distributed with the dispersion Uv, whereas the gas 
density is distributed log-normally with the mean no and the second central di¬ 
mensionless moment Uy (y = nn/uo)- Both stochastic fields are calculated using 
the Markovian processes (see LAK for mathematical basics). The set of param¬ 
eters estimated in the least-squares minimization of the objective function (see 
eqs.[29] and [30] in LAK) includes and Uy along with the total hydrogen col¬ 
umn density the mean ionization parameter Uq, and the metal abundances 
Za for a elements observed in a given absorption-line system. 

The computations are carried out in two steps; firstly a point in the pa¬ 
rameter space {A^h, Do, cJv, CTy, Za} is chosen at random and then an optimal 
configuration of {vi} and {rii} for this parameter set is searched for. These 
steps are repeated till a minimum of the objective function (y^ ~ 1 per de¬ 
gree of freedom) is achieved. To optimize the configurations of {vi} and {rii}, 
the simulated annealing algorithm with Tsallis acceptance rule (Xiang et ah, 
1997) and an adaptive annealing temperature choice is used (details are given 
in LACM). 


3. Main results 

From the above mentioned QSO spectra, we have selected 22 absorption systems 
suitable for the inversion analysis. Our results show that the QSO intervening 
systems are a heterogeneous population which is formed by at least three groups 
of absorbers : (1) extended metal-poor gas halos of distant galaxies; (2) gas 
in dwarf galaxies; and (3) metal-enriched gas arising from the inner galactic 
regions and condensing into the clouds within the hot galactic halo (galactic 
fountain). While the interpretation of a single system is sometimes subject to 
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Figure 1. {left panel) Carbon abundances plotted against the log¬ 
arithmic linear size of the corresponding systems as measured by the 
MCI procedure. Small size QSO absorbers seem to have systematically 
higher metal content than systems with larger sizes. 

{right panel) Plot of the line of sight velocity dispersion log((Tv) vs. 
log(A^H L) for the same sample of the metal systems. The dashed line 
is the linear regression log((Tv) oc k log(iVH L) calculated for the points 
shown as filled circles having both horizontal error bars. The slope n 
is equal to 0.30 ± 0.03 (1 a c.h). Open squares represent HVC-like 
systems. 


large uncertainties, the existence of a wide spread of properties in the different 
systems is certainly proved. 

Fj^ure 1 {left panel) shows an example of the measured carbon abundances 
[C/H]EI versus logarithmic sizes of the studied systems. Systematically higher 
metal abundances are seen in compact systems with linear sizes L ^ 0.4 kpc. 
The highest metallicity of [C/H] = 0.59 ± 0.08 was measured in the associated 
system at Zabs = 1.697 toward HE 0515-4144, the linear size of this system is 
only about 1 pc (Agafonova et ah 2003a). The lowest carbon content ([C/H] = 
—2.93 ±0.13) was found in the LLS at Zabs = 2.917 toward HE 0940-1050, which 
shows L ~ 140 kpc (Levshakov et al. 2003b). In this LLS we also measured 
[Si/C] = 0.35 ±0.15 and [Al/Si] = —0.18 ±0.13. These values certainly rule out 
the enrichment by very massive stars explosions and constrain the stellar masses 
by 25Mq. If, nevertheless, the Pop HI stars were massive (M > I 4 OM 0 ), 
then the average metagalactic metallicity produced by their explosions should 
be much lower than 10’ -3^©. 

If metal systems are formed in gas clouds gravitationally bound with inter¬ 
vening galaxies, the internal kinematics of the QSO absorbers should be closely 
related to the total masses of the host galaxies. In case of galactic population, 
different types of galaxies show different scaling relations between the linear size 


^Using the customary definition [X/H] = log(X/H) - logfX/H)©. Photospheric solar abun¬ 
dances are taken from Grevesse & Sauval (1998) and Holweger (2001). 
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and the velocity width of emission lines (e.g., Mallen-Ornelas et al. 1999). Pos¬ 
sible correlation between the absorber linear size L and its line-of-sight velocity 
dispersion was also mentioned in LACM. 

The correlation between and L stems from the virial theorem which 
states : ^ M/L no = N^L. Assuming that the gas systems are in 

quasi-equilibrium, one can expect ~ (iVHL)V2. 

In Figure 1 {right panel) we examine our systems by comparing their kine¬ 
matics (fJv) with measured sizes (L) and total gas column densities (A^h)- It 
is seen that in the log((Tv) versus log(A^H A) diagram, most systems with linear 
sizes L > 1 kpc he along the line with the slope k = 0.30 ± 0.03 (1 a c.h). 

Taking into account that we know neither the impact parameters nor the 
halo density distributions, this result can be considered as a quite good fit to the 
expected relation for the virialized systems. Hence we may conclude that most 
absorbers with L > 1 kpc are gravitationally bound with systems that appear to 
be in virial equilibrium at the cosmic time when the corresponding Lya absorbers 
were formed. If the most metal absorbers identified in the QSO spectra arise in 
the galactic systems then the question to what degree the intergalactic matter 
is metal enriched remains open. 

The absorption systems shown in Figure 1 {right panel) by open squares 
exhibit characteristics very similar to that observed for different types of HVCs in 
the Milky Way and may be interpreted as the high-redshift counterparts of these 
Galactic objects. The most clear example of such high-z HVC is the mentioned 
above very metal-rich system at Zabs = 1-697 with the neutral hydrogen column 
density A^(H i) = 4.4 x 10^^ cm“^. The system shows absorption from highly 
ionized transitions of C iii, C iv, N v, O vi. Si iv, and probably S vi. We found 
that only a power law ionizing spectrum {J^, oc is consistent with the 

observed sample of the line profiles, i.e., the system was probably blown out from 
the QSO/host galaxy by a cumulative effect of the supernova in the starburst 
and the quasar activity. Such metal-enriched outflowing gas is supposed to give 
rise to some HVCs observed in the Milky Way halo (Bregman 1980; Wakker 
2001 ). 

Some of the small-size metal systems show very strong metallicity gradients. 
A system at Zabs = 1-385 toward HE 0515-4414 is a good example. Here we 
measured V(H i) ~ 5.3 x 10^^ cm“^, [C/H] ~ —0.5, and L ~ 2 kpc. This system 
is probably embedded in an extremely metal-poor halo with V(H i) ~ 1.1 x 10^^ 
cm“^, [C/H] < —4, and L ~ 90 kpc (Agafonova et al. 2003a). The velocity shift 
between these subsystems is only Av = 140 km s“^. The velocity excess of the 
Galactic HVCs is usually greater than 90 km s“^, which is consistent with this 
Av value. 

In conclusion, we would like to emphasize that in spite of a utterly inho¬ 
mogeneity in metal abundances revealed in the Lya absorbers, all of them are 
perfectly described within pure photoionization models. Our results show that 
the fraction of the supposed shock-heated hot gas with temperature T > 10^ K 
is negligible in the analyzed absorbers, and that it is impossible to relate unam¬ 
biguously highly ionized absorption systems to the ‘warm-hot’ gas predicted in 
cosmological models (e.g., Cen &: Ostriker 1999; Dave et al. 2001). Each highly 
ionized system requires a comprehensive study in order to determine its nature 
and origin. 
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